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Aharonov-Bohm (AB) oscillations are studied for a parallel coupled vertical double quantum dot 
with a common source and drain electrode. We observe AB oscillations of current via a one-electron 
bonding state as the ground state and an anti-bonding state as the excited state. As the center gate 
voltage becomes more negative, the oscillation period is clearly halved for both the bonding and 
antibonding states, and the phase changes by half a period for the antibonding state. This result 
can be explained by a calculation that takes account of the indirect interdot coupling via the two 
electrodes. 

PACS numbers: 73.63.Kv, 73.23.Hk 



Quantum coherence and the correlation of electrons in 
semiconductor nanostructures are of primary interest. In 
the past few years, they have been well characterized for 
qubit states in single and double quantum dots (DQDs) 
in the context of quantum information. Before this, the 
direct observation of quantum coherence in QDs had al- 
ready been achieved using the Aharonov-Bohm (AB) ef- 
fect Conductance or current through an AB ring 
oscillates periodically with magnetic flux penetrating the 
ring area. The inclusion of a QD in the AB ring causes 
a phase change in the AB oscillations corresponding to 
changes in the number of electrons in the dot 0]- The 
abrupt phase change by the half period was initially ob- 
served [31 , and explained usingthe reciprocal theorem for 
a two-terminal measurement |3|. 

AB oscillations in current flowing through two QDs or 
DQDs embedded in an AB ring with one QD in each 
arm were studied theoretically with particular attention 
to the electron correlation effect The AB phase of 

a current flowing through two-electron states in a tunnel 
coupled DQD is predicted to be different by half a pe- 
riod in the cotunneling regime between the spin singlet 
state and triplet state due to interdot exchange coupling 
0. In contrast, for one-electron states, it is anticipated 
that the AB oscillations of the bonding state (BS) are 
phase changed by half a period from those of the anti- 
bonding state (AS), reflecting the parity of the electron 
wavefunctions in the DQD {q}^ . 

These features have been experimentally studied for 
parallel coupled lateral DQDs [10[ and vertical DQDs 11 1 
with a common source and drain electrodes, but without 
reference to the phase evolution. The phase change in the 
AB oscillation has only recently been observed in associ- 
ation with inelastic cotunneling processes but unrelated 
to interdot coupling [12|. In all previous reports, the ef- 
fect of interdot coupling was too weak to be resolved in 
the AB oscillations. 

In this paper, we use a laterally coupled vertical DQD 



with a common source and drain electrodes to study 
the AB oscillations of current through the one-electron 
states. We first precisely identify the one-electron BS as 
the ground state and the AS as the excited state in the 
charging diagram, and then measure the AB oscillations 
of current through the BS and AS. We observe critical 
changes in the phases and periods of the AB oscillations, 
depending on the center gate voltages. To account for 
these results, we consider the effect of the extra interdot 
tunnel coupling through the electrodes ("indirect inter- 
dot tunnel coupling"). 

Our device consists of laterally coupled vertical DQDs 
with four split gates made from a double-barrier het- 
erostructure (DBH) [HO]. A current flows verti- 
cally through the two parallel coupled QDs as shown in 
Fig. Ha). 

It is usually assumed that electron coherence readily 
collapses in the electrodes. However, this was not the 
case in our previous experiment on the same type of DQD 
device as used here because AB oscillations were indeed 
observed for the AB ring encircled through the source 
and drain ll|. This is because the electrodes are adja- 



cent to the thin AlGaAs tunnel barriers just outside the 
QDs. Therefore, we assume for the vertical QDs that the 
two QDs are coupled via the n-AlGaAs contact layers 
in a sufficiently short path length to maintain the coher- 
ence. Here, we introduce the indirect interdot tunnel cou- 
pling parameter and ad {\ois(d) \ = l)j which indicates 
the degree of coherence in the source and drain, respec- 
tively [8|. When as{d) = 0, the coherence is completely 
lost; the AB oscillation is not observed. In contrast, the 
coherence in the electrodes is completely conserved at 
|Q^s(d)l = 1; the visibility of the AB oscillation is at its 
maximum. Note that the visibility of the AB oscillation 
is 100% when the energy level in each dot is aligned. In 
actual systems, however, such a |Q:s(ji)| = 1 case is very 
special and most experimental situations correspond to 
|Q!s(d)l < 1- In particular, in the measured device, we 
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FIG. 1: (a) Diagram of current flow (Isd) in parallel cou- 
pled double quantum dot. The arrows indicate the current 
flow, (b) Color plot of the current as a function of the left 
and right side gate voltages, Vsl and Vsr, at a source drain 
voltage VsD=50 /iV and a center gate voltage Vc = —1.20 V. 
(Nl , Nr) denotes the electron numbers in the left and right 
dots, (c) Color plots of the average values of the transcon- 
ductance {dlso/dVsL 4- dIsn/dVsR)/\/2 as a function of Vsl 
and VaR at Vc = —1.20 V and Vsi5=800 ^iV, corresponding 
to the region within the circle in (b). 



assume |Q!s(d)| ^ 0.2, as described below [13 1 



Then two current loops with more or less the same 
area, S/2, are formed due to the presence of the direct 
interdot tunnel coupling. The magnetic flux, which pen- 
etrates the large current loop of the area S, including 
the two small same-sized loops, is defined by $ = BS as 
shown in Fig. [ija), where B is an externally applied in- 
plane magnetic field. Thus, the magnetic flux, which pen- 
etrates each small loop, is $/2. The evaluated AB oscil- 
lation period is ABq = 2$o/'S' = 0.77T, where $o = h/e 
is the magnetic flux quantum, h is the Planck constant 
and e the electron charge 13]. All the transport mea- 
surements are performed on a device placed in a dilution 
refrigerator with a base temperature of 20 mK. 

Figure Wljo) shows current (Isd) peaks or Coulomb 
peaks evolving with left and right gate voltages, Vsl and 
VsR, measured at a bias voltage, Vs£)=50 ^V. The cen- 
ter gate voltage, Vc, is fixed at -1.20 V. No Coulomb 
peaks are observed for Vsl ^ —1.7 V and Vsr £ —0.1 
V, indicating that the DQD is empty. Then, the elec- 
tron number -/Vl(A^_r) in the left (right) QD increases by 
one whenever the point in the parameter space in Fig. [T] 
(b) crosses one of the Coulomb peaks on the right (top) . 
Therefore, the charge states with (N^jN^) are fixed in 
each hexagon formed by the Coulomb peaks. 

Figure [IJc) shows the average value of the transcon- 
ductance {dlso/dVsL + dlso / dVsR) / ^ as a function of 
VsL and VsR at Vc = -1.20 V and 1/sd =800 ^V. The two 
anti-crossing stripes indicate that the Coulomb peaks are 
widened by applying a finite Vsd- The lower left, and up- 
per right stripes correspond to the one-, and two-electron 
states, respectively. The red or black lower border line 
of each stripe indicates the ground state, and the red or 
black line inside the lower left stripe is the one-electron 
excited state For the one-electron stripe, we observe 
that the ground and excited states anti-cross with each 
other as a typical feature of the resonance between the 
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FIG. 2; Current flowing through the bonding state (a), (c), 
(e) and (g) and the anti-bonding state (b), (d), (f) and (h) 
as a function of B for 14 =-1.20, -1.26, -1.36 and -1.40 V 
and at Vsd =800, 700, 500 and 500 /^V, respectively. Each 
inset shows the corresponding fast Fourier transform (FFT) 
spectrum of the current oscillation, where ABq = 2$o/S' = 
0.77 T. 



BS and AS, which are indicated by the circle and square 
in Fig. [TJc), respectively We estimated an interdot 
tunnel coupling 2t of 0.53 meV. The interdot Coulomb 
coupling energy Vmter is approximately 0.4 meV. The 
tunneling rates from the source and drain electrodes to 
the DQD, Ts and F^ are estimated to be ~ 0.5 /leV and 
~ 1.5 /ieV, respectively. 

For the two-electron state, the excited state is not well 
resolved in Fig. [T](c) , probably because the exchange cou- 
pling energy that corresponds to the energy separation 
between the singlet and triplet states is very small [l6j . 
Therefore, the state indicated by the triangle is the two- 
electron ground state. 

We measured the current at the circle and square 
points as a function of B for various Vc values. The two 
points are defined by the charging diagram measurement 
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at B =0T (not shown). Here, let us consider the BS and 
AS currents. The current at the circle in Fig. [Ijc) flows 
solely through the BS because only the BS is within the 
bias window, which is caused by applying Vsd- Con- 
sequently, the current is the BS current. On the other 
hand, at the square the current flows through the BS and 
AS, since both states are in the bias window. Using the 
BS and AS as the bases, the diagonal elements of a trans- 
mission matrix indicate the transmission probabilities of 
the BS and AS, and the total current is calculated by in- 
tegrating the transmission probabilities in the bias win- 
dow with respect to energy [H, 17[. Moreover, the cur- 
rent flowing through the system generally accompanies 
the mixing component of the BS and AS currents, which 
corresponds to the higher order tunneling processes com- 
prising both BS and AS. However, these processes, which 
are proportional to the second order terms of ajj^j), can 
be neglected, since as(d) is very small in our case |13|. 
Therefore, we can simply describe the total current as 
the sum of the BS and AS currents. Then, we can as- 
sume that the BS currents at the circle and square are 
almost the same, because the potential barrier height, 
i.e. the electron transmission probability, at the circle 
and square hardly changes. Consequently, we subtracted 
the current measured for the BS from that at the square 
to extract the AS current. 

The currents obtained for the BS, and AS are shown as 
a function of B in Fig. [2][a) , and (b) , respectively, where 
Vsd and Vc are the same as those in Fig. [UJc). Both cur- 
rents oscillate periodically with B; these oscillations are 
AB oscillations. The fast Fourier transform (FFT) spec- 
tra of (a) and (b) in the insets show peaks at 0.81 T, and 
0.79 T, respectively. These values are consistent with the 
calculated value of A_Bo=0.77 T. Therefore, this oscilla- 
tion corresponds to area S/2 rather than area S owing 
to direct interdot tunnel coupling 0, Q • There is a good 
half period contrast between (a) and (b) because there 
is a current dip at = T for BS whereas there is a 
peak for AS. It should be noted that the small dip at 
around ±1 T in Fig. [2] (b) is probably because the con- 
tribution of the BS current is not completely subtracted 
in the derivation of the AS current. 

The BS and AS currents obtained for Vc =-1.26, -1.36 
and -1.40 V and at Vsd =700, 500 and 500 /zV are shown 
in Fig. [JJc) to (h), respectively. Here 2t is estimated to 
be 0.46, 0.35 and 0.24 meV at Vc =-1.26, -1.36, and -1.40 
V, respectively. The interdot Coulomb coupling energy 
Vinter IS approximately 0.4 meV for the three different Vc 
values. In Fig.[2jc) and (d) at Vc = —1.26 V, the current 
oscillations are not very periodic but on the whole they 
are similar to those in Fig. ^a) and (b), respectively. 
The FFT spectrum in Fig. [2][d) shows that the largest 
peak is slightly different from that in Fig.ISfb), reflecting 
the larger peak spacing at S = T and ±1.1 T. 

The currents obtained for the more negative Vc values 
show markedly different features. At Vc = —1.36 V, the 



BS current oscillates only weakly with B in Fig. [UJe). 
The oscillation period of 0.5 T as identified in the FFT 
spectrum is approximately half the value obtained for 
the BS current at Vc = -1.20 V (a) and -1.26 V (c). In 
contrast, for the AS current in Fig.[5Jf), we observe clear 
oscillations with a period of 0.44 T, which is almost the 
same as that for the BS current in (e). The AB oscillation 
of the AS current in (f) with Vc — —1.36 V has a dip at 
B = Q T, whereas it has a peak at Vc — —1.20 V (b) 
and Vc = —1.26 V (d). On the other hand, that of the 
BS current in (e) for Vc = —1.36 V seems to leave a 
small peak at S = T, judging from the magnified plot. 
However, because the BS current is smaller than 1 pA for 
< 1 T and comparable to the noise level, it is difficult 
to resolve the AB oscillation at around B = T. 

The AB oscillations are significantly damped for both 
BS and AS, when Vc becomes even more negative at - 
1.40 V in Fig. M^g) and (h). The features of the AB 
oscillations appear similar to those at Vc — —1.36 V, 
although the AB oscillation of the BS is not very clear. 
The application of a more negative voltage to the center 
gate voltage weakens the direct interdot tunnel coupling. 
Moreover, it increases the dephasing in the electrodes, 
because of the increase in the effective path length in the 
electrode to connect the two dots. The latter could be 
the reason for the reduced visibility of the AB oscillations 
at Vc = -1.36 V and -1.40 V in Figs.lHe) to (h). 

However, possible changes in the distance between the 
two dots and the direct interdot tunnel coupling and in 
the dephasing caused by applying a much more negative 
Vc cannot account for either the abrupt changes in the 
phase or the period of AB oscillations with Vc as observed 
in Fig. [21 Here, we consider that the interference effect 
induced by the indirect interdot tunnel coupling via the 
electrodes [3] accounts for the observed features. 

First, let us simply analyze the currents flowing 
through the BS and AS, when the bias voltage is so large 
that the BS and AS are in the bias window. We de- 
fine the tunneling rate, F^^^^^ through the BS and Taa^\ 
through the AS from the source (drain) electrode to the 
DQD. They arc described as Q 
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rfr^ = r.,){i- 



■ as{d) cos(7r$/$o)} 
as(d) cos(7r$/$o)} 



(1) 
(2) 



where r^j^;-) corresponds to the tunneling rate from the 
source (drain) electrode to the DQD at as = ocd = 0, oth- 
erwise the Q!s(rf) parameters are assumed to be the same 
between the BS and AS. Because the higher order tun- 
neling processes comprising both the BS and AS can be 
neglected due to the small as(d) in our case as explained 
above, the currents, lb through the BS, and la through 
the AS are expressed as 

4(a) il/T^b{aa) + ^ f^bbiaa))^^ ^ (3) 

respectively [l3| . When as and have the same sign, 
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FIG. 3: Currents flowing through the bonding (solid line) and 
anti-bonding (dashed line) states as a function of magnetic 
flux $ at (a) as = —0.1 and ad = —0.2, (b) as = —0.1 and 
ad = 0.2 and (c) as — —0.1 and Od = 0.15, where F = Fs+Fd. 



lb {la) directly reflects the oscillatory parts in F^^ and 
^bb i^aa 8'iid r^jj) with the same period of 2$o in phase. 
Then, as expected from Eqs. ^ and ([2]), the AB oscilla- 
tion phase is different by a half period, i.e. $0i between 
lb and la- On the other hand, when as and ad have op- 
posite signs, the oscillation amplitude of lb (la) becomes 
small because F^j^ and F^ (Ff^ and F£,) osciUate with 
the same period of 2$o but in the opposite phase. Con- 
sequently, the AB oscillations of and la appear to be 
in phase with the period of 2$o but with an extra dip or 
peak, although the visibility of the AB oscillations should 
be lower than that when as and ad have the same sign. 
Therefore, it is possible to account for the changes in the 
period and phase of the AB oscillation with respect to Vc 
by using Us^d}- 

Here, let us discuss the experimental result in Fig. [2] 
concretely, using the change of Q;s(d)- Initially, we con- 
sider the signs of as and ad at Vc = —1.20 V in our 
experiment. The phases of the AB oscillations for the 
BS and the AS in Fig. ^a.) and (b) are not consistent 
with those predicted for as = = 1 This discrep- 
ancy can be explained by assuming that both as and 
ad have negative signs. Then, we consider the changes 
in as and ad as Vc becomes more negative. as((j) de- 
creases to with oscillation as the path length in the 
source (drain) electrode increases [13|. ad probably de- 
pends on the change in Vc more strongly than as, be- 
cause the screening effect is stronger in the drain than 
in the source in the side of the substrate. Therefore, for 
simplicity we assume that as remains a negative con- 
stant as tts = —0.1, and that ad changes from nega- 
tive to positive, i.e. ad = —0.2 to 0.2 then to 0.15 
and we calculate the currents flowing through the BS 
and AS numerically using the non-equilibrium Green's 
function method. The results are shown in Fig. [3] (a) 
to (c). We used 2t/hT = -265 and eVso/hV = 400 
in (a), 2t/hT = -175 and eVso/fiT = 250 in (b), 
and 2t/hr = -120 and eVso/hV = 250 in (c), with 
Vinter/fi^ ~ 200, kBT/KT = 10 as common parame- 
ters, where ks is the Boltzmann constant, T is the tem- 
perature, h = h/2TT and F = Fs + F^ (Fs = 0.3F and 
Yd = 0.7F), respectively. We evaluated all these param- 
eters experimentally. 



Figure[3ja) shows AB oscillations with the same period 
of 2$o but misaligned by $o between the BS and AS, so 
that the BS and AS currents have a dip, and a peak at 
$ = 0, respectively. These features are consistent with 
those at Vc = -1.20 V in Fig. dja) and (b). When ad 
becomes 0.2, which corresponds to the more negative Vc 
[l3| . Fig. [3l[b) shows almost the same AB oscillations 
with a period $o for both BS and AS; the period of the 
AB oscillation changes. Furthermore, the AS current at 
$ = in (b) is a dip, while it is a peak in Fig. Elja); 
the phase of the AB oscillation of the AS current also 
changes. These features are comparable to those in Fig. [2] 
(e) and (f) {Vc = -1.36 V). When ad becomes 0.15, 
which corresponds to a more negative 14, Fig. |3](c) also 
shows AB oscillations with a period 2$o for both BS and 
AS but with an extra small dip. This feature is similar 
to that for the AS current in Fig. [2] (h) at V;=-1.40 V, 
although the AB oscillation for the BS current is not clear 
in Fig. [2jg). Note that the features of the calculated AB 
osciUation change smoothly from Fig. [3l[a) to (c) (isf . 
Regarding Fig.[2ljc) and (d), because Vc — —1.26 V is an 
intermediate point between Vc = —1.20 and -1.36 V, ad 
probably represents the change from negative to positive. 
Therefore, the AB oscillations of the BS and AS are not 
as clear as those at Vc = —1.20 V because of the small 
absolute value of ad- Note that the small dip in the AS 
current at <& = in Fig. [Ijd) probably evolves into the 
large dip in Fig. [TJf). 

The AB oscillations of the BS are not as clear as those 
of the AS in Fig. [2Ue)-(h). We do not fully understand 
the reason for this, but consider that the Y s{d) value of 
the BS is smaller than that of the AS at around $ 
when as(d) is negative (see Eq.([T|) and ([U), resulting in 
the smaller BS current at around zero magnetic field. 

In conclusion, we have investigated AB oscillation for 
the BS and AS in a DQD. The periods of the AB os- 
cillations are roughly halved and, for the AS, the AB 
phase shifts by $0i depending on the center gate volt- 
age. These features are assigned to contributions from 
indirect interdot tunnel coupling, which depends on the 
center gate voltage. 
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DEVICE STRUCTURE 




FIG. SI: Schematic of double quantum dot device. The red arrow shows the direction of the apphed magnetic field B. 



Our device consists of two laterally coupled vertical quantum dots (QDs) with four split gates made from a double- 
barrier heterostructure (DBH) of an undoped 10 nm GaAs well and two undoped AlGaAs barriers with thicknesses 
of 8.0 and 7.0 nm as shown in Fig. Sl(a)) [l], [1]. Two of the four gates (side gates) are used to tune the electron 
number in each QD independently, and the remaining two gates (center gates) are used to modulate the strength of 
the interdot tunnel coupling. A current flows vertically through the two parallel coupled QDs between the common 
source and drain as shown in Fig. 1(b). 

The calculated area, 5/2, of the two current loops formed by the double QD (DQD) and either the source or drain 
electrodes due to the presence of the direct interdot tunnel coupling is approximately 350 nm x 15 nm= 5.3 x lO^^'^ 
m^, using our estimates of ^350 nm for the dot diameter and ~15 nm for the half value of the length of the double 
barrier structure. 
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INDIRECT COHERENT COUPLING PARAMETER "a" 
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FIG. S2: (a) Model of source (drain) electrode and two dots, s^j^) correspond to the distance between the two dots in the 
source (drain) electrode. (b)Dependence of a^^a) on length between two dots Ss(d). 

The indirect coupling parameter a{\a\ ^ 1) represents the degree of coherence in the electrodes [3|. Figure [S2] 
shows the calculation of as[d) as a function of the path length, Ss{d), in the source (drain) electrode to connect the 
two dots, when we simply assume one-dimensional electron propagation over the length Ss(d) as shown in Fig. IS2r b) 
0- kp is the Fermi wave number. As shown in Fig. IS2r a). as(d) decreases to with an oscillation as Ss(_d) increases. 
Note that the Ss(^d) dependence of Q!s(d) in (a) is stronger than that assumed in the calculation in Fig. 3. This is 
probably because the model (b) used for the calculation is too simple to provide a quantitative estimation of ag(^d) in 
the experiment. 



DERIVATION OF TUNNELING CURRENT THROUGH A DQD 

We derive the tunneling current through the DQD. For clarity of discussion, we concentrate on the noninteracting 
situation, and neglect the spin degree of freedom. In the last part of this section, the effect of the interdot Coulomb 
interaction is taken into account numerically. 

Using the nonequilibrium Green's function method, we obtain the following expression for the tunneling current 
through a DQD 



Here fu{e) is the Fermi-Dirac distribution function of the reservoir defined by 

Me) = ^ 



(1) 



(2) 



where fii, is the electrochemical potential of the reservoir v, and r(e, (f>) is the transmission probability defined as 

T(e,0) ^ Tr{G'-(6,0)r^((/))G''(6,0)r^(0)}, (3) 
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where the boldface notation indicates a 2 x 2 matrix whose basis is the tunnel-coupled bonding (6) and antibonding (a) 
states in a DQD. G'^^°'\e, ((>) is the retarded (advanced) Green's function matrix and r''(0) is the linewidth function 
matrix for the reservoir ly G S, D given by 



S /l + a^cos(7r$/$o) sin(7r$/$o) A /.x 

^"P' ~ ' \ -iasSm{n^/^o) 1 - a, cos(7r$/$o) / ^ ' 

-pD,,N P / 1 + arfcos(7r$/$o) -iad sin(7r$/$o) \ 
'^^^ mrfsin(^$/$o) l-adcos(7r$/$o) y ' 

r(</.)-r^(<^) + r^ (</.). (7) 

Here eb{a) is the energy of the tunnel-coupled bonding and antibonding states, and is the coherent indirect coupling 
for the reservoir u and (j) is the AB phase defined as = 27r<i>/<i>o, where $ is the magnetic flux threading through 
an AB interferometer and = h/e is the magnetic flux quantum. Moreover, we introduce the notation My , which 
denotes the matrix element of 2 x 2 matrix M. 

Using the above results, we obtain the following expressions for the tunneling current 

IsD{(t>) = h{(t>) + Ia{(t>) + hai<P), (8) 

hi^) = ^ 1 ^[/s(e)-/D(e)]rM|G^,(e,</))P, (9) 

lai^) = ^/ f [/5(e)-/D(6)]rirfjGL(e,0)P, (10) 
hai^) = ^1 f [/5(e)-/D(e)][rf,rfjG:;,(6,0)p-frf„rf,|GL(e,0)P 

+2Re {rf,rf,G^,(e, 0)[GS,(e, 0)]* + rf,rf,G^,(6, ^)[GL(e, </>)]* 
+rf,rf,GL(e, 0)Kj,(e, 0)]* + rf,rf,G^,(e, 0)[GL(e, 0)]* 

+rirS,GL(e,0)[GL(e)]* + rf,rf,G;,(6,0)[GL(e,0)]*}] • (n) 

Here we evaluate /&(</>) in a high bias voltage regime such as eVso ^ — — "^tc, ^ij{4')ikBT 

hW = [I + Oia^a,,)] , (12) 

J- bb 

where tc is the direct interdot tunnel coupling. In our paper, we focus on weak coherent indirect coupling conditions, 
namely |ai,| 1. We assume that the terms proportional to a^a^> are negligible. Then, we have 



t bb 
1 



1 4- 1 



(13) 



This is equivalent to Eq. (3) in the main text. This result has the form of a Breit-Wigner type transmission through 
the bonding state. Therefore, we call Ibi4') the bonding state component of the tunneling current. Similarly, we can 
derive the expression of the antibonding state component of the tunneling current. 

Ia{ct>) ^ e , (14) 

Moreover, the mixing component of the bonding and antibonding states for the tunneling current is given as 

76a(0) oc 0{a„a„>). (15) 
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FIG. S3: Numerical results for the tunneling current, (a) Bonding and antibonding components, (b) Mixing component. 



Thus, this mixing component is negligible when we use the above approximation. 

The above results are the same even when there is a finite interdot Coulomb interaction. Under the condition shown 
in Fig. 3 (b) in the main text, we provide numerical results for the bonding, antibonding, and mixing components 
for the tunneling current in Fig. IS3I It is clear that the mixing component Iba{4') of the tunneling current is much 
smaller than the bonding Ib{4') and antibonding /a('/') components of the tunneling current. Moreover, we find that 
the mixing component vanishes when '^/'^q is an integer. This can be explained as follows. According to Eqs. ([5]) 
and (|ni), the off-diagonal matrix elements of the linewidth functions are zero when 4>/$o is an integer. In this case, 
the bonding state is orthogonal to the antibonding state. Thus, the bonding and antibonding states do not contribute 
mixing component, namely Iba{4') = 0- 



as AND ad DEPENDENCES OF AB OSCILLATION OF BONDING AND ANTIBONDING STATES 



We calculate the bonding and antibonding state currents at an intermediate ad value in Fig. lS4r a) between Fig. 3(a) 
and (b), and in Fig. [S4F b) between (b) and(c), respectively. We used 2t/hT = —200 and eVso/HT = 250 in (a) and 
2t/hT = -150 and eVso/fir = 250 in (b), respectively, with Vinter/KT = 200 and ksT/hT = 10 as the common 
parameters. 

In Fig. 3(a), the period of the osciUations for the bonding and antibonding states is 2$o with a dip at $ = 2ri$o 
($ — (2n + l)$o) for the bonding (antibonding) state, where n is an integer. As ad increases to ~ 0.14, the oscillation 
for the bonding (antibonding) state has an additional small dip at $ = {2n+ l)<I>o = 2n$o) as shown in Fig. lS4r a). 
Thus, the period of the AB oscillations for the bonding and antibonding states changes to $o but with two alternating 
amplitudes. The small dip grows and becomes comparable to the other dip when ad ~ 0.2 as shown in Fig. 3(b). As 
ad decreases to ^ 0.18, the difference between the two alternating oscillation amplitudes increases again (see Fig. IS4I 
(b)), and finally, the AB oscillations for the bonding and antibonding states evolve into those seen in Fig. 3(c). 



AHARONOV-BOHM OSCILLATION FOR TWO-ELECTRON GROUND STATE 



The coherence, i.e. the AB oscillation, of the two-electron state is attracting interest to detect the singlet and 
triplet states similarly to the that of the one-electron states. Consequently, we investigate the AB oscillation 
for the two-electron ground state. Figure ISSl a) and (b) show experimental data obtained at Vc =-1.20 and -1.36 V, 
respectively. We observe periodic current oscillations with a large peak at i? = T in (a) but no clear oscillations in 
(b). The FFT spectrum in (a) indicates two oscillation periods of 0.81 T and 1.8 T. For the first period, we find the 
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FIG. S4: Currents through the bonding (soUd Une) and antibonding (dashed line) states as a function of a magnetic flux $, 
for the indirect coupling parameter in the source and drain electrode, (a) cvs = —0.1 and aa = 0.14 and (b) = —0.1 and 
aa = 0.18 




FIG. S5: Current flowing through the two-electron state at (a) Vc =-1.20 V and (b) Vc =-1.36 V as a function of B. The insets 
show the FFT spectra of the current oscillations. 



same period in Fig. 2(a) and (b) but not for the other, which arises from the three large peaks at and ±1.8 T in 
Fig.EHa). 

We used a master equation method expanded from that reported in [6| to calculate the current flowing through the 
two-electron ground state. 

When the exchange coupling energy is very small the singlet and triplet states are not well resolved. Then, the 
Heitler-London (HL) state comprising equivalent contributions from the bonding and anti-bonding states is the two- 
electron ground state. The current through the HL state does not oscillate with the magnetic field because the 
oscillations of the bonding and antibonding states are canceled out [6j. Therefore, we can obtain the AB oscillation 
in the relatively strong tunnel coupling case as shown in Fig lSST a). 

On the other hand, when the exchange coupling is strong enough for the singlet and triplet states to be well 
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separated, the singlet state is composed of the HL and molecular states for which two electrons with anti-parallel 
spins occupy the bonding state. Therefore, the current through the molecular state oscillates. Accordingly, we cannot 
obtain the AB oscillation in the weak tunnel coupling case as shown in FiglSSfb). 

The two phases of the AB oscillations in Figs. [S5] (a) and (b) are the same as those of the anti-bonding state as 
shown in Figs. 2(b) and(f), respectively. However, this is not consistent with the theory, because the AB oscillation 
phase of the molecular state is the same as that of the bonding state 0. This discrepancy can be understood if as 
and ad for the two-electron state are different from those for the one-electron state. This argument is reasonable, 
because the Vsl and Vsr values of the one-electron and two-electron states are different, even when Vc is the same. 
A more detailed study is necessary for the two electron ground state. 
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